A Unified treatment of small and large- scale dynamos in helical turbulence 
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Helical turbulence is thought to provide the key to the generation of large-scale magnetic fields. 
Turbulence also generically leads to rapidly growing small-scale magnetic fields correlated on the 
turbulence scales. These two processes are usually studied separately. We give here a unified 
treatment of both processes, in the case of random fields, incorporating also a simple model non- 
linear drift. In the process we uncover an interesting plausible saturated state of the small-scale 
dynamo and a novel analogy between quantum mechanical (QM) tunneling and the generation of 
large scale fields. The steady state problem of the combined small/large scale dynamo, is mapped 
to a zero-energy, QM potential problem; but a potential which, for non-zero mean helicity, allows 
tunnehng of bound states. A field generated by the small-scale dynamo, can 'tunnel' to produce 
large-scale correlations, which in steady state, correspond to a force-free 'mean' field. 



PACS Numbers : 47.65. +a, 52.30.-q, 95.30. Qd 

Large-scale magnetic fields in astronomical objects are 
thought to be generated by dynamo action involving heli- 
cal turbulence and rotational shear ||l|,^ . Here large-scale 
refers to scales much larger than the outer scale, say L, 
of the turbulence. However, turbulent motions, with a 
large enough magnetic Reynolds number (MRN hence- 
forth), can also excite a small-scale dynamo, which ex- 
ponentiates fields correlated on the tubulent eddy scale, 
at a rate much faster than the mean field growth rate 
These two dynamo problems, viz. the small-scale 
dynamo (SSD) and large-scale dynamo (LSD), are usu- 
ally treated separately. However this separation is often 
artificial; there is no abrupt transition from the field cor- 
related on scales smaller than L and that correlated on 
larger scales. We show here that the equations for the 
magnetic correlations, which involve both the longitudi- 
nal and helical parts, are already sufficiently general to 
incorporate both small and large-scale dynamos in the 
case of random fields. They provide us with a paradigm 
to study the dynamics in a unified fashion, which could 
be particularly useful to study possible inverse cascade 
of magnetic fields, to scales larger than L. 

Consider the induction equation for the magnetic field. 

(aB/9t) V X (v X B-77V X B), (1) 

where B is the magnetic field, v the velocity of the fluid 
and 77 the ohmic resistivity. Take v = vt + v^, the 
sum of an externally prescribed stochastic field v-r, and 
a drift component V£), which models the non-linear back 
reaction of the growing Lorentz force. We assume to 
be an isotropic, homogeneous, Gaussian random velocity 
field with zero mean. For simplicity, we also assume vt 
to have a delta function correlation in time (Markovian 
approximation). Its two point correlation is specified as 
< i;^(x,i)4(y, s) >= T'^{r)6{t - s) with (cf. §) 

r^(r) = T^[S^^ - (^)] + nC-^) + Ce^.frf. (2) 



Here <> denotes ensemble averaging, r = |x — y| and 
= — y*. TLir) and T^ir) are the longitudinal 
and transverse correlation functions for the velocity field 
while C (r) represents the helical part of the velocity cor- 
rellations. If V.vt = 0, Tat = {l/2r)d{r^TL)/dr. (Note 
vt is also not correlated with the magnetic field.) 

To model the drift velocity in a tractable manner, but 
one which nevertheless gives some feel for possible nonlin- 
ear effects, we proceed as follows. As the magnetic field 
grows, the Lorentz force pushes on the fluid. We assume 
the fluid almost instantaneously responds to this push, 
and develops an extra, 'drift' component to the velocity 
proportional to the instantaneous Lorentz force. So we 
take a model v^i — a[(V x B) x B)], with the parameter 
a = r/(47r/9), where r is some response time, and p is 
the fluid density. (Such a velocity can also arise when 
friction dominates inertial forces on ions, as in ambipo- 
lar drift). This gives a model nonlinear problem, where 
the nonlinear effects of the Lorentz force are taken into 
account as simple modification of the velocity field. Such 
a phenemenological modification of the velocity field has 
in fact been used by Pouquet et al. and Zeldovich et 
al. (pg. 183) [|l| to discuss nonlinear modifications to the 
alpha effect, and we adopt it below. 

Consider a system, whose size S >> L, and for which 
the mean field averaged over any scale is zero. Of course, 
the concept of a large-scale field still makes sense, as the 
correlations between field components separated at scales 
r » L, can in principle be non-zero. We take B to be a 
homogeneous, isotropic, Gaussian random field with zero 
mean. This is a perfectly valid assumption to make in 
the kinematic regime (v^i = 0), as the stochastic vt, 
has these symmetries. When we include the non-linear 
drift velocity, it amounts to making a closure hypothesis, 
which we do here again for analytical tractability. The 
equal-time, two point correlation of the magnetic field is 
given by < B^{x,t)B^ (y^t) >= AP^{r,t), where 
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= Mn[5'' - (^)] + Ml{^) + He 



(3) 



(Here <> denotes a double ensemble average, over both 
the stochastic velocity and stochastic B fields) . Ml {r, t) 
and Mn{t, t) are the longitudinal and transverse correla- 
tion functions for the magnetic field while H{r^ t) repre- 
sents the (current) helical part of the correlations. Since 
V.B = 0, Mat = {l/2r)d{r'^ML)/{dr). 

The stochastic Eq. (|l]) can be converted into the evolu- 
tion equations for and H. We give a detailed deriva- 
tion of these equations elsewhere, including the effect of 
the non- linear drift ||]. We get 



dt dr ^ '^^ dr ' 
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where we have defined 

kn^V + Tl (0) - Tl (r) + 2aML (0, t) 
UN = 2C(0) - 2C(r) - 4ai/(0, t) 
G = -4 [(TAT/r)' + {rTLy/r'] 
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Here prime denotes a derivative with respect to r. These 
equations form a closed set of nonlinear partial differen- 
tial equations for the evolution of Ml and H, describ- 
ing, as we will see, both SSD and LSD action for ran- 
dom fields. The effective diffusion kn includes micro- 
scopic diffusion (rj), a scale-dependent turbulent diffu- 
sion {Tl(0) — Tl^t)) and nonlinear drift adds an amount 
2aMi(0, t), proportional to the energy density in the fluc- 
tuating fields. Similarly aN represents a scale dependent 
a-effect (2C(0) — 2C(r)) and nonlinear drift decreases 
this by 4a_ff(0, i), proportional to the mean current he- 
licity of the magnetic fluctuations. This modification to 
the a-effect is the same as that obtained in [^]^. The 
G{r) term allows for the rapid generation of magnetic 
fluctuations by velocity shear and the existence of a SSD 
independent of any large-scale field . 

First consider non-helical turbulence, with C{r) ~ 0, 
allowing solutions with H{r,t) = 0. This case has been 
extensively studied for the kinematic case (a = 0), when 
TL{r) has a single scale (cf. ||l|,p| ) and by us for a model 
Kolmogorov type turbulence g7 These studies show that 
one has SSD action, and magnetic fields correlated on 
scales upto the turbulent scale can be generated. In the 
kinematic limit, kjv and G are time independent. One 
can then look for eigenmode solutions to (jj), of the form 
*(r)exp(2rt) = t'^^/k^Ml- This transforms Eq. (|) 
for ML{r,t), into a time independent, Schrodinger-type 
equation, but with a variable (and positive) mass, 



d^^ 

r* = -KA,— -|-C/o(r)«'. 



The "potential" is Uo{r) = T'l + {2T'Jr) -h k^/2 - 
(K'Ar)^/(4KAr) -|- 2knIt'^ , for a divergence free velocity 
field. The boundary condition is '5 — > 0, as r ^ 0, oo. 
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Note that 

TL{0))/r'^ a.s r 00. The possibility of growing modes 
with r > obtains, if one can have a potential well, with 
Uq sufficiently negative in some range of r, to allow the 
existence of bound states, with an "energy" E = —T < 0. 

Suppose we have turbulent motions on a single scale 
L, with a velocity scale v. Define the magnetic Reynolds 
number (MRN) Rm = vL/t]. Then one finds that 
there is a critical MRN, = i?c ~ 60, so that for 
Rm > Rc, the potential Uq allows the existence of bound 
states. For R^ = R^, F = 0, and this marginal sta- 
tionary state is the "zero" energy eigenstate in the po- 
tential Uq For Rm > i?c, r > modes of the 
SSD can be excited, and the fluctuating field correlated 
on a scale L, grows exponentially, on the corresponding 
'eddy' turn-over time scale, with a growth rate Fi ~ v/ L 
0. To understand the spatial structure of the flelds, 
define w{r,t) B(x, t).B(y, t) >= {l/r^)d[r^ML]/dr, 
the correlated dot product of the random field. For the 
SSD, w(r) is strongly peaked within a region r — rd ~ 
^Rm)"^/"^ about the origin, for all the modes, and for 
the fastest growing mode, changes sign accross r ^ L and 
rapidly decays with increasing r/L. (For the marginal 
mode with F = 0, Rm is replaced by Rc). Note that rd 
is the diffusive scale satisfying the condition rj/r'^ ^ v/L. 
A pictorial intepretation of the correlation function, due 
to the Zeldovich school (cf. 0,^), is to think of the field 
as being concentrated in "flux ropes" with thickness of 
order rd << L, and curved on a scale up to ~ L, to 
account for negative values of w. 

How does the SSD saturate? The back reaction, in 
the form of a non- linear drift, simply replaces rj by an 
effective, time dependent tjl) = t] + 2aMi(0,i), in the 
Kj\i term of Eq. (jj). Suppose we deflne an effective 
MRN, for fluid motion on scale L, by Roit) = vL/rjD{t). 
Then as the energy density in the fluctuating fleld, say 
Esit) = 3Afi(0, t)/87r, increases, Rl> decreases. In 
the flnal saturated state, with (dML/dt) — (obtain- 
ing say at time t^), Ml, and hence the effective rjo ui 
(^ become independent of time. Solving for this sta- 
tionary state then becomes identical to solving for the 
marginal (stationary) mode of the kinematic problem, 
except that Rm is replaced by RL>{ts)- The flnal sat- 
urated state is then the marginal eigenmode which ob- 
tains, when Eb has grown (and Rd decreased) such that 
Rnits) = vL/{r] + 2aMi(0,ts)) ^ Rc 60. Also w{r) 
for the saturated state will be strongly peaked within a 
region r = Vd ~ L[Rc)^^^^ about the origin, change sign 
accross r ^ L and then rapidly decay for larger r/L. 
From the above constraint, AlL{0,ts) — vL/{2aRc), as- 
suming T] << 2aML{0,ts). So at saturation. 
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Note that r is an unknown model parameter. If we were 
to adopt T ~ L/v, that is the eddy turn-over time, then 
Eb at saturation is a smaU fraction ~ i?^^ ^ 1, of the 
equipartition value. Further suppose we intepret w{r) 
for the saturated state in terms of the Zeldovich et al. 
picture; of flux ropes of thickness r^, curved on scale L, 
in which a field of strength Bp is concentrated. In this 
picture, the average energy density in the field Eb ^ 
{B^/8Tr)LryL^. Using r^/L^ « R-^, and r - L/v, we 
then have Bp/Sir ^ /3w^/2, where, remarkably, the R^^ 
dependence has disappeared, and Bp has equipartition 
value. So the SSD could saturate with the small-scale 
field having peak values of order the equipartition field, 

— 1/2 

being concentrated into flux ropes of thickness LRc , 
curved on scale L, and an average energy density R~^ 
times smaller than equipartition. 

We now turn to consider the effect of helical correla- 
tions. If ao = 2C(0) 7^ 0, then one can see from Eq. 

and that new generation terms arise at r >> L, 
due to the a- effect, in the form Ml — ■■■■ — AarH and 
H = ... + [r^iaTMLyy/r"^. Here ar ^ ao - AaH{0,t) 
and dot represents a time derivative. These couple Ml 
and H and lead to the growth of large-scale correla- 
tions. There is also decay of the correlations at r ^ L, 
due to diffusion with an effective diffusion co-efficient, 
rjT = Tj + Tl{0) + 2aMi(0,t). From dimensional anal- 
ysis, the effective growth rate is Tl> ~ ar/D — rir/D^, 
for correlations on scale ~ D, as in the large-scale a^- 
dynamo. This also picks out a special scale Dq ^ rjT/fXT 
for a stationary state (see below). Further, as the SSD, 
is simulataneously leading to a growth of Ml at r < L, 
in general at a faster rate v/L >> ao/D, the growth of 
large-scale correlations can be seeded by the tail of the 
SSD eigenfunction at r > L. The SSD generated small- 
scale field can thus seed the large-scale dynamo. Indeed, 
as advertised, both the SSD and LSD operate simula- 
taneously when aQ ^ 0, and can be studied simply by 
solving for one ML{r, t). 

The coupled time evolution of H and Ml for a non-zero 
Q!o requires numerical solution. But interesting analytical 
insight into the system can be obtained for the marginal, 
stationary mode, with [dML/dt) = (dH/dt) = 0. Infact 
both the kinematic and non-linear dynamo problem can 
be treated in a unified fashion. With H independent 
of time, Eq. (||) implies 2knH + uj^iMl = 0, for any 
solution regular at r = and vainishing at r — s- cx) . From 
this, as r ^ 0, 2r/i?(0, t) = 0, and hence (0, i) = for a 
non-zero ?7, a result which also follows directly from the 
evolution (conservation) of magnetic helicity H]. So any 
general nonlinear addition to the a-effect which arises in 
terms of H{0,t) has to vanish in a stationary state! 

Now substitute H{r) = -aAr(r)Mi,(r)/(2KAr(r)) into 




r > L 



FIG. 1. Schematic illustration of the potential U (r) for the 
marginal mode in helical turbulence. A non zero uq allows 
the tunneling of the zero-energy state to produce large-scale 
correlations. 
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and define once again 



fK^ML- We get 



We see that the problem of determining the magnetic 
field correlations, for the marginal/stationary mode once 
again becomes the problem of determining the zero- 
energy eigen-state in a modified potential, U ^ Uq — 
4[C(0) - C(r)]2/Kjv. Note that the addition to Uq, due 
to the helical correlations, is always negative definite. So 
helical correlations tend to make bound states easier to 
obtain. When C(0) = 0, and there is no net a-effect, 
the addition to Uq vanishes at r » L, and U 2riT/r^ 
at large r, as before. The critical MRN for the station- 
ary state, will however be smaller than when C(r) = 0, 
because of the negative definite addition to Uq. 

When QfQ — 2C(0) 7^ 0, a remarkable change occurs in 
the potential. At r 3> -/j, where the turbulence velocity 
correlations vanish, we have U{r) — 2riT/r'^ — a^/rjT- 
So the potential U tends to a negative definite con- 
stant value of —a^/rjT at large r (and the effective mass, 
1/2kn l/2rjT, a constant with r.) There are strictly 
no bound states, with zero energy/growth rate, for which 
the correlations vanish at infinity. We have schematically 
illustrated the resulting potential U in Figure 1, which is 
a modification of figure 8.4 of Zeldovich et al. 0| . In fact, 
for a non-zero a^, U corresponds to a potential which al- 
lows tunneling (of the bound state) in the corresponding 
quantum mechanical (QM) problem. It implies that the 
correlations are necessarily non-zero at large r > L. The 
analytical solution to (|^) at large r 3> is easily ob- 
tained. We have for r L, ML{r) — ML{r), 

^^^('') = ^372 [Ci-^3/2(A'^) + C2 J_3/2(H], (10) 

where fj, = ao/r/T = Dq^ , and Ci, C2 arbitrary constants. 
Also w{r) = w{r) — ^r^^[Ci sin/ir -I- C2 cos/ir]. Clearly 
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for a non-zero ao , the correlations in steady state at large 
r, are like "free-particle" states, extending to infinity! 

An alternate derivation of ML{r), for the kinematic 
case, clarifies its meaning further. Suppose one thinks 
of the large-scale field as a "mean" field Bq, the mean 
taken over cells much larger than L. Assume Bq itself 
is random over different cells, statistically homogeneous 
and isotropic, with a correlation < BoiBoj >= M^{r). 
Let M^, and be respectively the corresponding 
longitudinal, transverse and helical correlations. Then 
Bq in each cell obeys the kinematic, mean-field dynamo 
equation (dBo/dt) = V x (aoBo - (?/ -I- Tl(0))V x Bq), 
whose steady state solution is, V x Bq = ^ioBo, where 



imposed on Mlj, 



jj,o = ii{a = 0). This constraint 
gives = and (Mf 

filML/2. Also as V.Bq = 0, = {l/2r){r^MLy. 

These three equations fix all the functions uniquely. We 
get, remarkably, M^{r) = ML{r), with /i replaced by the 
kinematic value /xq- So this solution actually describes a 
random mean-field, for the marginal large-scale dynamo. 
Similarly if we had imposed V x Bq — /iBo, would be 
given by Eq.(|l^). Note this also shows that the effective, 
steady-state, large-scale field Bq is force-free, although 
B itself is not. 

It is straight-forward to connect the large-scale, force- 
free field for the marginal mode of helical turbulence, 
with the SSD generated field, as they are both the so- 
lution of the same Eq. (^), for large and small r re- 
spectively. For example, one can integrate Eq. (||), 
adopting different starting values of Ml{0, tg), and taking 
M']^{0,ts) = 0, to construct a whole family of solutions 
(parameterised by Mi(0,ts)), which match small-scale 
correlations with the large-scale correlation of (|o|). For 
each such solution, we will have one value of C1/C2. Note 
that this is unlike the standard QM tunneling problem, 
where the boundary condition that the free-particle state 
is an outgoing wave at large r uniquely fixes the tunnel- 
ing amplitude, for a given bound state. However, when 
we consider a zero-energy, stationary state, there is no 
such natural time-asymmetric boundary condition; so no 
unique fix for C1/C2, in ®. Nevertheless, if AfL(0,i,) 
is so small, and Rd so large, that U admits bound states 
with energy E < —o^/tit, then the corresponding time 
dependent system is unlikely to lead to stationary cor- 
relations. This sets a lower bound on Mi(0,ts) or Eb- 

Further as ML(0,t) = + lQH'^{Q,t), the saturation 

of ML(0,i), depends on how fast H{0,t) 0, its sta- 
tionary value; over and above the saturation effects of 
the increasing diffusion rjjy. So the full time dependent 
problem needs to be solved to fix an upper bound on 
Mi(0,ts) or Eb |10|. We will return to this elsewhere. 

In conclusion, we have given here a unified treatment 
of small- and large-scale dynamos, for the case of random 
fields, incorporating also a simple model non-linear drift. 
We uncovered an interesting plausible saturated state of 



the small-scale dynamo. For random fields, we argued 
that any non- linear addition to the a effect in terms of the 
average current helicity H{0, t) (cf. has to vanish in 

a steady state. The steady state problem of the combined 
small/large scale dynamo, was then mapped to a zero- 
energy, QM potential problem; but a potential which, 
for ao ^ 0, allows tunneling of bound states. A field 
generated by the SSD, can then "tunnel" to seed the 
growth of large-scale correlations, which in steady state, 
correspond to a force- free "mean" field. It remains to 
solve the fully time-dependent system and to incorporate 
more realistic back-reaction effects of the Lorentz force. 
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